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ABSTRACT 

ERATO is a Monte Carlo event generator describing four fermion production at 
LEP2 and beyond. All tree-order processes leading to four fermions are included, 
taking into account all relevant QCD contributions. QED higher order corrections 
are introduced via a structure function approach, whereas weak corrections related 
to fermion loops are also included. Special attention has been paid in studying the 
trilinear gauge couplings, and all possible contributions, including the CP-violating 
ones, have been implemented. 
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PROGRAM SUMMARY 



Title of the program: ERATO. The name originates from Sgarcv, the ancient Greek Muse of 
lyric poetry. 

Program obtainable from: Dr. Costas G. Papadopoulos, Institute of Nuclear Physics, NRCPS 



'Democritos', 15310 Athens, Greece and using ftp: // alice.nrcps.ariadne-t.gr / pub / papadopo/ erato 
Licensing provisions: none 

Computer for which the program is designed and others on which it has been tested: HP, IBM, 
ALPHA and SUN workstations. 

Operating system under which the program has been tested: UNIX 
Programming language: FORTRAN 77 

Keywords: four-fermion processes, trilinear gauge couplings, event generator 

Nature of physical problem: A very important fraction of the events produced at high-energy 
e + e~ collisions corresponds to four-fermion final states. Important physical issues, like the mea- 
surement of the mass of the W boson and the study of the trilinear gauge couplings (TGC), are 
based on the analysis of these four-fermion final states. Other production processes, like asso- 
ciated Higgs production or R-parity violating SUSY particle production, lead also to the same 
final states. It is therefore indispensable to have a rather accurate description of the four-fermion 
production, including tree-order signal and background contributions as well as the leading part 
of the higher order corrections. 

Method of solution: The construction of an event generator is of course the desired solution 
to the problem of calculating all four-fermion processes at high energies. To this end, we have 
calculated all relevant matrix elements, using a variation of the spinor technique which is more 
efficient in writing and testing the corresponding FORTRAN codes. A multichannel approach 
on phase space generation has been used in order to deal with the problem of the multipeak 
structure of the integrated function. Higher order corrections are systematically introduced and 
the new physics effects described effectively by the trilinear gauge couplings are considered. 
ERATO provides the computational framework which enables us to study in a systematic way 
four-fermion physics at high-energy e + e~ colliders. 
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LONG WRITE-UP 



1 Introduction 

LEP2 is already functioning at the highest energy ever reached for e + e~ colliders. Moreover, 
machines capable to produce even higher energies, up to the TeV scale, like a Next Linear 
Collider (NLC), are nowadays thoroughly studied Most of the physicists working in the 
field think that new physics is to be expected at these energies and a lot of work has been 
done more than two decades now trying to develop theoretical frameworks to incorporate the 
expected 'new phenomena'. On the other hand our current understanding of the elementary 
particle interactions, the so called Standard Model , is not yet completely tested due to the 
undetected Higgs particle a very important ingredient of the present theoretical apparatus. 
Any new dynamics in the Higgs sector of the theory will have very important consequences 
on the interactions among the heavy vector bosons, W and Z |2|, H|. In order to study these 
consequences, an accurate knowledge of the heavy vector boson production at these energies is 
indispensable and since heavy vector bosons always decay into light fermions, the four-fermion 
production is undoubtly the first most interesting channel(s) to be studied at these energies. On 
the other hand, four-fermion will be a background channel for most of the interesting expected 
'new phenomena', like SUSY particle production, new heavy vector boson contributions, etc. 

ERATO is an event generator, which enables one to have an accurate description of the 
four-fermion dynamics. In this paper we will try to give a presentation of this program, starting 
with the main theoretical features and ending with a short description of how it works with a 
few illustrative examples. In section 2 a brief presentation of the amplitude and phase space 
calculations is given, in section 3 the main features of the FORTRAN code are described and finally 
in section 4 certain calculations are presented. 

2 Computational Framework 
2.1 Helicity Amplitudes 

The first important step in order to construct an event generator is of course the calculation of 
all relevant tree-order matrix elements contributing to the processes under consideration. Our 
first approximation consists of neglecting the fermion masses, for all external fermion legs: this 
is of course a rather good approximation, since these terms are of the order of mj/s, where mj 
is the mass of the light fermion and s is the square of the center of mass energy. Nevertheless, 
this reasoning fails in the case of small angle scattering, a rather special case for which an 
approximation will be described in the sequel. In ERATO we have used a special representation 
of the so called 'spinorial structure' of the amplitudes, which follows more closely the structure 
of the Feynman graph under consideration, thus facilitating the writing and testing of the code. 
More specifically we are using the so called E- vector formalism^, ||, which is based on the 
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following representation of the fermion 'current', for massless fermions: 
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with p^ = p° ± p 3 




Figure 1: Feynman graphs contributing to e + e — > e u e qq' . 

Let us see how the Feynman graph of Fig.||(a), contributing to the process e~ (pi)e + (P2) — * 
^(^3)^(^4)9(^5)9^ '{p&) 1 is computed in this formulation. First of all the spinorial part is simply 
written as: 

u~{p- ) )^{-p , G-h-p , A)l Vu ~{P&) u-ipa^u-ipi) u x {p2)l^u x (pi) 
= - b i E -(P5,Pi) • E-(P3,P<0 E_(pi,p 6 ) ■ E x (p 2 ,Pi) (3) 

i=3,4,6 

where 6, = ±1 depending on whether the particle is outgoing or incoming and A = ± corresponds 
to the helicity of the incoming electron. Then coupling constants and propagators are combined 
to give the final expression. The graph of Figjl|(b) is simply computed by the interchange of 
P2 *-* P3 and 62 <-> 63. All garphs have been computed using the above framework. Moreover 
where graphs involving as external particles massless bosons, as in the case of e~e + — ► qqgg 
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which contributes to the four-jet production channel and is experimentally indistinguishable 
from four-quark production, the polarization vector of the outgoing gluon is given by: 

e»(p;\)^^=EZ(p, q ) (4) 

The auxiliary four-vector g M reflects the existence of the gauge invariance, always accompanying 
massless spin-1 particles, and provides a very powerful test of the computation, since any choice 
°f Q (q " P ¥^ 0) should give exactly the same result. Further details on the computation of the 
relevant color matrix can be found in references || 0]. 

2.2 Phase Space Generation 

The next important part is of course the phase space generation. This follows closely the 
reasoning of references []|, [|. The problem is that the amplitude we have to integrate over is 
a very complicated function of the kinematical variables, peaking at different regions of phase 
space. The idea is to define different kinematical mappings, corresponding to different peaking 
structures of the amplitude and then use an optimization procedure to adjust the percentage 
of the generated phase-space points, according to any specific mapping, in such a way that the 
total error is minimized. Let us present a very simple example in order to clarify these issues. 
The differential cross section is written as: 

da = ^-\M\ 2 dUps{P; Pl ,..., PA ) (5) 
As 

where 

n n 

c£Lips(P;pi, ..., Pn ) = (2n) 4 ~ 3n ]J (Vp* S(pf) 0(p°)) S(P - J» . (6) 

i=l i=l 

For the double resonant graph of Fig.§ and in order to describe efficiently its peaking structure 
we have naturally chosen the variables s + = p\ = (p q + Pq>) 2 and s_ = p 2 _ = (pi + Pu e ) 2 and 
write the invariant phase space as 

dlA^P^p^pu^Pq^Pq') = dLips(P;p+,p_) dLips(p + ;pq,Pq>) ^ dLips(p-;pe,pp e ) ■ (7) 

In this way we have now the freedom to generate the integration variables s± such that their 
distribution follows the well known Breit-Wigner form, namely 

1 



(s± - m 2 w f + m 2 w Y 2 ¥ 
This is actually achieved in ERATO by the following mapping: 

m \q< = m w + rnw^w tan(pi (y + - y_ ) + y_ ) 

, ~i(E 2 — rn w . 

V+ = tan ( 

mwi- w 

V- = tan (-— ) 
1 w 

m tu, = m w + rnw^w tan(p 2 (y+ -y-)+y~) 
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Figure 2: Double resonant graph contributing to e + e 



y' + 

cos Ow 
cos Q\ 
cos6* 



m 



2/95 — 1) 4>£ = 2vr/96 (in the £, 9^ rest frame) 
2p7 — 1, 0* = 27rps (in the g, g' rest frame) 



(8) 



where pi 



, 8, represent uniformly distributed pseudorandom numbers in (0,1). In order 



to have a rather complete description of the process e + e~ — * e~v e qq', for example, sixteen 
different mappings have to be used. The number of mappings used is specific to each process 
and is appropriately defined in ERATO. 

The large number of mappings, whose phase-space densities are described by the probability 
distributions gi(£l), naturally introduces the following representation of the total phase-space 
density: 

g(n) = Y,Ui9m (9) 



where f2 stands for the phase space element and 

J dn g(tt) = J dn gi (Q) = a* = 1 



(10) 



The variables CKj, are known as the a priori weights. Their main property is that the result of 
the MC integration is independent of them: 



JV 



dnf(n) = (w)= iim 



(11) 



On the other hand the MC error depends on af. this enables us to perform a weight optimization 
which could increase the speed of the MC integration^]. This is achieved by minimizing the 

X A detailed description can be found in reference M. 
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variance estimator, W(a) = (w 2 ), with respect to ai, under the constraint X)« a i = 1- The 
optimization therefore is equivalent to the evaluation of the parameter values a° pt , for which 
the MC error is minimized. This leads to the following set of equations: 

Wi(a) = Wj{a) , for all i,j = l,...,N 

where 



with the MC weight w defined by 



Wi{a) = \m w ) 



ERATO uses an iterative algorithm to solve the above equations. The convergence of this 
procedure is checked by the evaluation of the variable 

V = max\Wi-WA (12) 

which should vanish at the optimum point. In the numerical procedure T> measures indeed how 
well our a's approximate the optimum solution ol^ 1 . 

2.3 Higher Order Corrections 

A substantial part of the higher order corrections to the tree-level amplitudes have been taken 
into account in ERATO. As is well known the main part of these corrections comes from the 
emission of soft photons from the initial colliding particles. Although for neutral current pro- 
cesses, this ISR correction is a well defined, gauge-invariant quantity, for charged currents only 
the leading logarithmic (LL) part can be computed unambiguously. Nevertheless, this LL part 
accounts for most of these corrections, and therefore still it can be used to evaluate the ISR ef- 
fect. In ERATO we have implemented the ISR in the structure function approach. This simply 
means that the true differential cross section is not given by Eq.(||), but by: 

^°"ISR = dxidx2f(xi)f(x 2 )da(xiX2s) (13) 

where the momenta of the initial particles are given by, p e - = xi^-(l; 0, 0, 1) and p e + = 
12^(1; 0, 0, —1). The explicit form of the function fix) fl(| is: 



- + x)- l -f3 2 



I _|_ 3 X 2 

log(x) + 4(1 + x) log(l — x) + 5 + x 

1 — x 



(14) 



where ■je is the Euler constant and 



6 



with a being the electromagnetic coupling constant. 

An other important part of the higher order corrections is the so-called Coulomb correction. 
This is rather important at the threshold region, y/s ~ 161 GeV , reaching the level of 5%. It is 
the result of the long-range character of the electromagnetic interactions, contributing a factor 
roughly proportional to a//3, which near threshold gives a substantial contribution. Its explicit 
form is 



with 



°Coul - Born og 



1 



2 \Pm + *\ 2 -P 
1 arctan = — — ; — 

7T 1 2 /5Im(/3 A f) 



(15) 



P = -\Js 2 -2 S (p 2 + p 2 _) + (p 2 + -p 2 _) 2 



Pm = \/l-4:M 2 /s, M 2 = m 2 w - im w T w - ie 

A = ^~ P -[ (16) 
s 

and — 7r/2 < arctan y < tt/2. Here (3 is the average velocity of the W bosons in their centre-of- 
mass system and CC3 refers to the three, double resonant graphs contributing to the on-shell W 
pair production, e~e + — > W~W + . 

Finally an important part of the radiative corrections is related to the width of the unstable 
particles, W and Z bosons. As it has been shown in references JH, [llj the corresponding cor- 
rections, which consist of resuming all closed fermionic-loop contributions to the two and three 
point functions, although necessary to restore gauge invariance of the calculation, is rather small 
at LEP2 energies, with the exception of final state topologies where an outgoing electron (or 
positron) is very close, essential parallel to the beam axis. 

In order to describe the small angle scattering of processes with an e~ or e + in the final state 
we have to take into account the small electron mass, m e , leading to a further complication of the 
computational procedure. Nevertheless a very good estimate of the total massive cross section 
can be obtained by a rather simple approximation. This is achieved, at the leading logarithmic 
level by introducing a cut on the angle of the outgoing, massless electron (or positron) given 
by: 

_m e (p° - q°) 



p°q° 



(17) 



where p° (q°) is the energy of the incoming (outgoing) electron. In all cases we have checked, 
this approximation gives a rather good estimate of the total massive cross section. 

We conclude this section by underlining the fact that ERATO has historically been designed 
for TGC studies and followed an older MC for single W production j| . Moreover in the present 
version of ERATO, all TGC deviations, including CP violating ones, are considered^. 

3 Program Structure 

The complete file system is organized as follows: 

2 For a detailed presentation of TGC interactions, see reference Jl2[ 
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flag IPRO 


processes 


type of graphs 


1 




NC6 


2 


v T v r v T v T 


NC12 


3 


V T V T \X + \x~ 


NC10 


4 


fl + jl~ T + T~ 


NC24 


5 


T + T - T + T - 


NC48 


6 




NC20 


7 


e + e~r + T~ 


NC48 


8 




NC144 



Table 1: Physical processes in program lllLn. 



flag IPRD 


processes 


type of graphs 


1 


T + T~ V T V r 


NC12+CC7 


2 


e + e~v e T> e 


NC24+CC14+CC18 


3 


T + T~V e U e 


NC12+CC9 


4 


v T v T v e v e 


NC6+CC6 


5 


v e v e v e v e 


NC18+CC18 



Table 2: Physical processes in program MLcl. 



1. f files contains the seven FORTRAN files (* . f ) where the computation of the corresponding 
processes is performed. More specifically, tables [l||7| give the physical processes which can 
be obtained by each FORTRAN program. Also shown is the corresponding value of the flag 
IPRO which is used as an input. When more than one physical processes are assigned to a 
given value of the IPRO, it means that all these processes share identical description. Finally 
the number of Feynman graphs contributing to the current procces is also presented. The 
nomenclature NC24, CC10, etc stands for Neutral or Charged Currents and is related to 
the unstable states contributing to a given process. 

2. inputs contains the corresponding input files used in the actual version. 

3. demo contains the corresponding Make file used for compilation and linking as well as a 
test run output. 

4. share contains all commonly used subroutines and functions, such as the function EVECTOR 
and the pseudo-random number generator (s). 
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flag IPRO 


processes 


type of graphs 


1 


fJb~D li T + U T 


CC9 


2 


e~v e T + i> T 
e+v e ^~v^ 
e + v e T~D T 


CC18 



Table 3: Physical processes in program llll_c2. 



flag IPRO 


processes 


type of graphs 


1 


v^v^DD 
v T v r DD 


NC10 


2 


v^v^UU 
v r v r UU 


NC10 


3 


fj, + fj,~UU 

T+T-UU 


NC24 


4 


H + n~DD 
t+t-DD 


NC24 


5 


e+e-UU 


NC48 


6 


e+e~DD 


NC48 



Table 4: Physical processes in program llqq. 



The codes are very flexible and easily accessible to the user, so that any update can easily 
be implemented. It should be noted however that no special FORTRAN-optimization has been 
performed in the writing of the code. Nevertheless the time needed for event generation is 
rather small, so for all usual applications such an optimization is not really necessary. 

Let us start with the main common variables used in the program. The first of course 
is P (1:4, 1:20), where all particle momenta are stored: P (4, 1:20) refers to the energy and 
P(l:3,l:20) to the three-momentum, p x , p y , p z . At the beginning of each program (e.g. 
llll_n.f), the momentum assignment of the final particles is explicitly given, e.g. 

(ps)^iJ,{p4)u(pg)d(pio). Variable P(l:4,l) refers to the incoming electron and P(l:4,2) 
to the incoming positron. Included in the same common MOM is the array B(l:20) which as 



flag IPRO 


processes 


type of graphs 


1 


v e v e UU 


NC12+CC7 


2 


v e v e DD 


NC12+CC7 



Table 5: Physical processes in program neneqq. 
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flag IPRD 


processes 


type of graphs 


1 


t~v t UD 


CC10 


2 


e-p e UD 


CC20 


3 


H+VpUD 
t + v t UD 


CC10 


4 


e+v e DU 


CC20 



Table 6: Physical processes in program evud. 



flag IPRD 


processes 


type of graphs 


1 


UUDD 


CC11+NC24+QCD8 


2 


UU'D'D 


ecu 


3 


UUUU 


NC48+QCD16 


4 


DDDD 


NC48+QCD16 


5 


UUU'U' 


NC24+QCD8 


6 


UUD'D' 


NC24+QCD8 


7 


DDD'D' 


NC24+QCD8 


8 


UUgg 
DDgg 


QCD8 



Table 7: Physical processes in program qqqq. 



explained in the previous section, takes the values ±1 if the particle is outgoing (incoming). 
Physical constants are included in the common PHYS. 

The main routine where the helicity amplitudes are calculated, is called MASTER. The ar- 
ray WTI (1 : 2 , 1 : NITER) gives the amplitude squared, for =p helicities of the incoming electron, 
whereas NITER refers to different input values for non-standard TGC, with ITER=1 always re- 
turning the Standard Model value. MASTER is called by the routine EVENT after the phase-space 
point generation has been performed. This is done by first producing the appropriate energy 
fractions XI and X2, when ISR is on (flag ISR=l), which define the reduced c.m.s. energy, 
s = X1X2S0. Notice that the generation of these variables follows the function FISR(X) which 
is the product of the corresponding structure functions given by Eq.(14). Then the routine 
ADDRESS is called, which returns the phase space point, according to the given phase-space map- 
ping, labeled by the flag IGEN. The flag IFLAG is used as follows: when IFLAG=0, the phase space 
point is generated whereas for IFLAG=1 the phase space density is computed. Of course, as 
explained in the previous section, the distribution among the different mappings is governed by 
the values of c^, ALPHA (1 : NGEN) , where NGEN is the total number of mappings used in the actual 
calculation. Optimization on the values of oti, which means redefinition of them, is performed 
after certain iterations defined by the variable IBASE. EVENT is finally returning the value of 
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the variable WTI (1 : 2 , 1 : NITER) , which corresponds to the MC weight of the given phase-space 
point (event). The driving routine is called DRIVE, which takes the weight coming from EVENT 
and then performing all relevant sums. A schematic representation of the flow of the program 
is given in Fig.0. 



DRIVE 



ISR 

Xl,X 2 



EVENT 



W: 



ADDRESS 



CCUTS 

any cut 



MASTER 
\M\ 2 



ALGOR's 
Qi -> ^ 



EVECTOR 

Ex(p,q) 



Figure 3: The flow-chart of the program. 

For the cases where WW^f and WW Z vertices appear, as for example in the evud.f code, 
the flag INAIVE is used to distinguish between different ways of treating the width of the unstable 



particles. As explained in detail in reference [11|, the fixed width solution, corresponding to the 
INAIVE=1 can be safely used by default. INAIVE=0 and INAIVE=2 correspond to the 'fudge-factor' 



schemes 13], whereas INAIVE=5 to the running width, without corrections to the WW"f and 
WW Z vertices, which of course at small momentum transfer gives inconsistent results. Finally 
INAIVE=3 correspond to the inclusion of all the imaginary-parts of the fermion-loop corrections, 
in a way consistent with both £7(1) and SU(2) gauge invariance. 
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In order to avoid singularities of the massless amplitude as well to be as close as possible 
to the experimental picture cuts are applied in the usual sense. In the present version we have 
implemented the so-called Canonical Cuts, used in the LEP2 workshop analyses. Of course any 
cut can be trivially implemented. In the case of ISR, there is an additional cut on the reduced 
energy of the e + e~ system, defined by the variable SCUT. 

All matrix elements have been tested against the MadGraph |l4j calculations to sixteen digits 
accuracy. Also consistency checks related to U{\) and SU{2) gauge invariance, or equivalently 
to the high-energy unitarity have been performed. Finally ERATO has been participated to 
the comparisons tests made at the LEP2 workshop [15|, where a very accurate check on the 
generator whole has been successfully performed. 



4 Test Run 



A typical input file is as follows: 

1 *process 

1 *iterations 

1 *ISR 

*IC0UL0MB 

128.07 0.2310309 91.1888 2.4974 80.23 2.033 *input parameters 

*cuts (total xs) 
1.0 1.0 *cmax(cmin) cmas 
175 *energy 

200000 *nev 

2500 *IEV 

1 *0 no optimization 

1 *naive 
1.0 2. 1.5 2.0 0.8 0.5 0.0 

0.0 1. 0.5 1.0 -0.2 -0.5 -1.0 

0.0 1. 0.5 1.0 -0.2 -0.5 -1.0 

2 *icase 



This is read by the evud program. The process under consideration is defined by the flag IPR0. 
The NITER corresponds to how many iterations the code will perform in order to calculate the 
weight for NITER different values of the TGC parameters (default=l). Then ISR defines the use 
of initial state radiation (on=l, off=0). The same for IC0UL0MB which is the case of the Coulomb 
correction. The input parameters correspond to l/a em , sin 2 9 W , mz, Tz, mw, respectively. 
The flag ITOTAL, in case it takes the value 1, will use the approximation described by Eq.(|l7]) for 
computing an estimate of the total massive cross section. CMAX(CMIN) and CMASS are variables 
used by the generator algorithms and can be put to their limits, CMAX=1 and CMASS=0 respectively. 
Then the energy and the number of MC iterations (NEV) are given. Variable IEV defines the first 
time when optimization will take place and the flag I0PT switch on (off) this procedure. INAIVE 
governs the width scheme to be used, with INAIVE=1 being the fixed width, and INAIVE=5 the 
naive running width. Several other options are available for specific processes, which can be 
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found in the distributed codes. Then NITER values for TGC parameters with Standard Model 
values 1 and respectively are read by the program (the remaining terms are not read). Finally 
the flag ICASE pick up the desired TGC parameter-relation scheme^, known as ob$ (ICASE=0), 
aw® (ICASE=2) and a\y (ICASE=3) as well as other choices which can easily be found in the 
distributed subroutine MASTER. 

The output, apart form monitoring printings, looks as follows: 

ISR= 1 
IC0UL0MB= 1 

STARTING AT Wed Jul 31 22:18:49 1996 

THIS IS THE STARTING POINT 
ITOTAL = 

CMAX = 1.0 CMASS = .5 
E= 175.0 
NEV= 250000 

1 . /ALPHA , ZL0 , ZETA , GEUL , GAB 

137.036 25.53137686693119 5 . 698195932780109E-02 .5772156649015319 

.9701624225753434 

E0,ELMAS,ZETA1,SCUT 

175.0 5.000000000000000E-04 17.54941409170019 1.0 
RM48 INITIALIZED: 

CANONICAL CUTS 1.0 3.0 .984807753012208 .9961946980917455 25.0 1 

ICASE= 2 

IMPR0VMENT FOUND 

DISTANCE= 6.855963070506428E-06 

EVENTS 2500.0 2500 

IMPR0VMENT FOUND 

DISTANCE= 1.797322661221867E-06 

EVENTS 7500.0 10000 

IMPR0VMENT FOUND 

DISTANCE= 1.019962512017920E-06 

EVENTS 10000.0 20000 

IMPR0VMENT FOUND 

DISTANCE= 7.980700337882255E-07 

EVENTS 20000.0 40000 

IMPROVMENT FOUND 

DISTANCE= 6.970462448481489E-07 

EVENTS 40000.0 80000 

IMPROVMENT FOUND 

DISTANCE= 6.011073515995931E-07 

3 For a detailed analysis on these TGC parameters, see reference |p.6[. 
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EVENTS 80000.0 160000 



ALPHA ( 


1) = 


.000 


ALPHA ( 


2) = 


.298 


ALPHA ( 


3) = 


.006 


ALPHA ( 


4) = 


.003 


ALPHA ( 


5) = 


.003 


ALPHA ( 


6) = 


.006 


ALPHA ( 


7) = 


.683 



ENERGY = 175.0000 
ANOMALOUS MAG. MOMENT OF W= 

.OOOOOE+00 1.0000 1.0000 .OOOOOE+00 .OOOOOE+OO 



SIGMA= .4875463E-03 +- . 1390471E-05NB 



SIGMA= .3478089E-05 +- . 1243032E-07NB 

TOTAL SIGMA ERROR 
.491024E-03 . 138875E-05 

ICASE= 2 

[TOTAL, PASS, FAIL, ICCUT] : 250000241945 8055 4789 1479 815 972 



GENERATOR 784 78252 3938 3030 2968 4178 156850 000000000000 


[MAXIMUM WT=] 1.489366913311956E-02 
ENDING AT Wed Jul 31 23:56:49 1996 

This is a run for the process e~e + — > i^v^ud. The variables SIGMA give the left and right 
handed helicity contributions with respect to the initial electron. The TOTAL SIGMA and ERROR 
are the total cross section and its MC error. The flag IMPROVMENT FOUND refers to the result of 
the optimization and the final values of Oj are also printed out. The variable DISTANCE refer to 
V as defined in Eq.(|l^). Finally the values of the canonical cuts, Ef, Ej et , cos0£, cosd^ j^ and 
m'je.t-jet are presented after the flag CANONICAL CUTS. All the other elements in the output are 
more or less self-explanatory. 
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